Sequential quantum cloning under real-life conditions 
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We consider a sequential implementation of the optimal quantum cloning machine of Gisin and 
Massar and propose optimization protocols for experimental realization of such a quantum doner 
subject to the real-life restrictions. We demonstrate how exploiting the matrix-product state (MPS) 
formalism and the ensuing variational optimization techniques reveals the intriguing algebraic struc- 
ture of the Gisin-Massar output of the cloning procedure and brings about significant improvements 
to the optimality of the sequential cloning prescription of Delgado et al [Phys. Rev. Lett. 98, 150502 
(2007)]. Our numerical results show that the orthodox paradigm of optimal quantum cloning can in 
practice be realized in a much more economical manner by utilizing a considerably lesser amount of 
informational and numerical resources than hitherto estimated. Instead of the previously predicted 
linear scaling of the required ancilla dimension D with the number of qubits n, our recipe allows a 
realization of such a sequential cloning setup with an experimentally manageable ancilla of dimen- 
sion at most D = 3 up to n — 15 qubits. We also address satisfactorily the possibility of providing 
an optimal range of sequential ancilla-qubit interactions for optimal cloning of arbitrary states under 
realistic experimental circumstances when only a restricted class of such bipartite interactions can 
be engineered in practice. 

PACS numbers: 03.67.Lx; 03.67.Bg; 03.65.Ta; 02.70.-c; 71.27. +a 



I. INTRODUCTION 

The superposition principle of quantum mechanics 
precludes quantum information from being perfectly 
"cloned" and transformed into classical information. Such 
a fundamental impossibility and chief difference between 
classical and quantum information is formally phrased 
in terms of the so-called "no-cloning" theorem [1| of 
quantum computation and quantum information Q. 
Nonetheless, quantum cloning machines (QCM) allow 
copying of arbitrary and a priori unknown quantum 
states with an imperfect albeit optimal fidelity as a mea- 
sure of the quality of copies (clones) p-Q- Such QCMs 
provide an optimal way of isometrically (unitarily) evolv- 
ing an initial set of N input qubits all in the same un- 
known state (augmented by M—N blank qubits) into the 
final state of M approximate clones, with the whole evo- 
lution described by the map N —> M (M — > M). How- 
ever, engineering such maps through a single application 
of a global isometry (unitary) operation that is capable 
of entangling all input qubits simultaneously is in general 
prohibitively difficult to realize in an experiment. As a 
consequence, all experimental realizations of QCMs hith- 
erto either in optical systems [6h"12| or NMR setups [iH - 
Il5| have suffered from the requirement for exceedingly 
high experimental demands that particularly poses ma- 
jor obstacles in scaling the devices to multiqubits. In an 
effort to work around such an experimental challenge, a 
sequential prescription for an experimentally manageable 
implementation of the multiqubit entangler of the opti- 
mal QCMs has been proposed by Delgado et al [16j. In 



sequential quantum cloning paradigm of Delagdo et al 
an ancillary system is introduced to interact sequentially 
and only once with each qubit in a row and is set to de- 
couple from the qubit chain in the last step [16]. This 
is essentially equivalent to a nontrivial decomposition of 
the multiqubit entangling operation of the QCM into a 
one-way sequence of consecutive two-body ancilla-qubit 
operations (see Fig. [I). Although a general genuinely 
entangling multiqubit operation cannot be implemented 
via such a sequential prescription [13 13; the case of se- 
quential implementation of an N — > M isometry required 
for QCM tasks was proved to be always possible |17| . 
From experimental point of view, various physical set- 



(a) 



U 



12... M 



M-N r- 
blankj - 
qubits] - 

I 

input < - 
qubits 



QCM 




U 



Ma 



Figure 1. M —¥ M quantum cloning machine implemented 
(a) globally and (b) sequentially with the aid of an ancillary 
system. Each bipartite unitary Uka in (b) acts only on the 
Hilbert space of qubit k and ancilla a and leaves other qubits 
intact. The fact has been illustrated by using vertical lines 
with open circles at each qubit they act upon. Ancilla degrees 
of freedom are shown by solid bold lines throughout. 
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tings employ sequential operations for generation of en- 
tangled multiqubit states. Paradigmatic setups include 
the photonic streams of cavity or circuit QED with a 
cavity mode as the ancillary system [l9l - [26| and the se- 
quential generation of electron-spin entangled states in 
quantum dot setups utilizing a completely mixed nuclear 
spin bath as the ancilla [231 - In terms of the practical 
benefits, sequentiality and the absence of measurements 
guaranteed by a unitarily decoupling ancilla may tremen- 
dously facilitate the physical implementation of a gener- 
ically complex global cloning operation. 

Furthermore, from theoretical perspective, it has been 
demonstrated [16| that the multipartite entangled state 
of the output clones produced by a sequentially imple- 
mented quantum doner can be characterized in terms 
of the hierarchy of the so-called matrix-product states 
(MPS) f28] that had previously arisen in the context of 
spin chain models [29| and are by now recognized as a uni- 
fying framework for rephrasing and optimizing various 
numerical renormalization group techniques of strongly 
correlated systems (30l - l32l |. An MPS representation of 
an n-qubit output of a QCM is then given by 

1 

|Vout>= £ (ipFlVfa.-.V^lin,...^), (1) 
i n ...ii—0 

where Dk x _Dfc + i-dimensional matrix VjS represents the 
physical interaction between ancilla and kth qubit with 
local space \ik), and with \ipj) and \<pf) denoting the 
initial and final ancilla state, respectively. The bond di- 
mension of an MPS is defined as D = max^D^. 

As pointed out by Delgado et al, it turns out that the 
minimal bond dimension of the MPS representation of 
([T]) coincides the required ancilla dimension D (e.g., the 
number of atomic levels) for sequential realization of the 
desired quantum doner |16]. Thus, apart from the algo- 
rithmic facilitation of the sequential prescription in terms 
of the decomposition of a global doner into consecutive 
two-body interactions, its practical feasibility and scala- 
bility yet hinges crucially upon ancilla demands and the 
way it scales with the number of clone qubits. A de- 
tailed analysis of such an issue revealed that the required 
minimal ancilla dimension D scales only linearly with the 
number of clones M (more precisely D = 2M for the case 
of universal symmetric cloning of Gisin-Massar as de- 
scribed in the subsequent section). This result promises 
an experimentally affordable ancilla-assisted scenario for 
quantum cloning of a small number of qubits, but yet 
infers the one that demands more and more experimen- 
tal resources upon increasing the number of qubits. On 
the other hand, cloning of arbitrary unknown states en- 
tails access to a vast variety of ancilla-qubit interactions 
(V^ in Eq. ([T])), whereas some of the required ones may 
be unattainable within a given physical setup. In this 
sense, an important experimental question will be if the 
existing physical setups could still fulfill the resource re- 
quirements of the proposed sequential QCM (SQCM). 



And if the answer is no, to what extent the numerical 
resources of the SQCM algorithm can be further opti- 
mized toward less costly experimental demands and, in 
turn, allowing realization of SQCM under such restricted 
experimental conditions. It is the purpose of the present 
paper to address these issues. 

This paper is organized as follows: Sec. |TT] sets the 
scene by rephrasing the optimal cloning machine of Gisin 
and Massar in terms of MPS and providing some tech- 
nical details on simulation of such a doner. In Sec. IIIII 
the issue of restrictions on the ancilla dimension D is 
addressed and numerical "regularization" techniques for 
achieving the cloning task with the minimal possible re- 
sources are presented. Section ITVl deals with restrictions 
on ancilla-qubit interactions and optimization protocols 
for realizing the cloning task with the minimal accessi- 
ble class of physical operations. Finally, Sec. fV\ contains 
our conclusions and an assessment of the applicability 
of our proposed methods to future realization of scalable 
quantum cloning. 



II. MATRIX-PRODUCT REPRESENTATION OF 
THE GISIN-MASSAR STATE 

We consider here the universal symmetric cloning 0, 0| 
of a single qubit in an arbitrary unknown state \tp) — 
a\0) + (3\1) to M clones described by the isometry map 
1 ->■ M 

\1>) -> |VWt) = \GM M (rl>)) 

M-l 

= £ aj\(M-j)^,ji; x )s 

3=0 

®\(M-j-i)r,jr x )s , (2) 

with the coefficients aj — \J m(m~-S ^)'> anc ^ — 
j)4>, j(f>' L )s denotes the normalized completely symmetric 
state with M—j qubits in state \<p) and j qubits in the or- 
thogonal state \4> ± ). The resulting multiqubit state from 
the cloning procedure, the so-called Gisin-Massar state, 
describes an entangled state of M clones augmented by 
M — 1 anticlones that are introduced to guarantee the 
optimality of the cloning procedure. We shall later ex- 
ploit such a "parity" feature of clones and anticlones for 
encoding aj into MPS representation of the Gisin-Massar 
state. 

An (2M — l)-qubit state of Gisin-Massar written in 
computational basis of the form 

\GM M m= J2 E •••E c L™- 2 ,..., 1 

«2M-1=0 l2M-2=0 il=0 

|*2M-Xi*2M-2) ■ ■ ■ ,h) ,(3) 

can be cast into a canonical MPS representation with 
minimal bond dimension D (2ll . [23| . The symmetriza- 
tion requirement of the Gisin-Massar state within either 
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Figure 2. (a) Graphical representation of the MPS descrip- 
tion of Gisin-Massar state and (b) the overlap of states used 
to calculate the contractions in the cost function Eq. I|12|l. 
The boxes represent Vj fc j matrices of a sequentially prepared 
state of Gisin-Massar form. The links connecting the boxes 
represent indices that are being contracted (or summed over) . 
The more information (e.g. larger ancilla dimension or larger 
number of ancilla-qubit couplings) the matrices contain, the 
darker and larger the associated boxes are drawn throughout. 



clone or anticlone subspace establishes the algebraic con- 
nection between otj and . It indeed makes most of the 
coefficients vanish or to be identical. Distinct values 
of Q,j correspond only to a priori symmetrized products 
of clone kets and the anticlone ones and are equal in 
number to the number of clones M to be produced. We 
shall see in the subsequent section that this feature pro- 
foundly affects the simulability of such a state in terms 
of the required numerical resources. 

The first step in achieving an MPS representation of 
the Gisin-Massar state is to perform Vidal decomposi- 
tion [33[ on ([3|) . Treating a bipartite decomposition of the 



coefficient C 



as a matrix with indices 



~i2M — l\i2M — 2).")il 

*2M-I an d the combined one %im-i = *2Af-2, ■ • ■ , U , we 
perform singular value decomposition (SVD) f3~Il. |351| on 
C 



c = usv^ 



(4) 



Identifying now the left unitary as IA = ^M-l]' we P ro ~ 
ceed with the singular value decomposition of the remain- 
ing part 7Wl 2M_1 l = SV^ within the next bipartite de- 
composition Z2M-1, «2M-2 \i2M-3, ■ ■ ■ ,h and obtain 



M [2M-1] =U ' S ' V '\ 



(5) 



Identifying again the left unitary as U' = V[2M-2] anc ^ ^' 
crating such a procedure till the end of the chain yields an 
open-boundary MPS representation of the Gisin-Massar 
state of the form ([1]) with n replaced by 2 M — 1. A 
useful graphical representation of Gisin-Massar state in 
MPS form has been shown in Fig. Ha). 

Some remarks on technical details of reshaping the 
(or equivalently ajs) into ^-matrices are in order. An 
operational ambiguity may arise from contributions de- 
generate in otj coefficient, though, belonging to cloning 
of |0) or |1) 



M-l 

\GM M (0))= aj\(M-j)0,jl) s 

3=0 

®\(M-j -1)1, jQ) s , (6a) 

M-l 

\GM M (l))= a 3 \(M-j)l,jO) s 

®\(M-j -1)0, jl) 8 , (6b) 



whose linear combinations constitute the expression 
Eq. |2) for cloning of an arbitrary state It will be, 
however, essential to distinguish such degenerate contri- 
butions in each bipartite combination of the indices in 
order to realize an in situ restructuring of the informa- 
tion. For this purpose, we suggest a parity-bit encoding 
scheme to get around this difficulty. The scheme relies on 
the basic observation that contributions associated with 
cloning of |0) (|1)) enjoy even (odd) parity in terms of 
the number of qubits in state |1). Hence, such a degen- 
eracy is formally lifted upon the action of a parity-bit 
operator V with eigenvalues ±1 for a register of qubits 
with even (odd) parity. The technical details of the lat- 
ter procedure go beyond the scope of the present paper 
and will be published separately. The parity bit encod- 
ing scheme has already been employed in the context of 
classical and quantum cryptography [36], redundant ar- 
ray of independent disks (RAIDs) and optical quantum 
computing schemes (37j . 

Following such a constructive recipe, results in a Gisin- 
Massar state in canonical form with the maximal bond 
dimension D = dL"/ 2 J w here d denotes the local dimen- 
sion (d = 2 for qubits) [28]. Such a maximal bond di- 
mension will be the minimal required one for the canon- 
ical description of the target state, too, provided that all 
the constituent V matrices are full-rank, otherwise rank- 
deficient matrices could give rise to a bond dimension 
lower than the canonical one. In this sense, as we shall 
see in the subsequent section, rephrasing Gisin-Massar 
in terms of MPS allows to obtain vivid insights into the 
resource requirements of the Gisin-Massar cloning ma- 
chine. 
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Figure 3. (Color online) Results for the regularization of the ancilla dimension D. The inset of (a) shows the error (1-J-) on 
a log scale which turn out to be of the order of the numerical noise (depicted by the dash-dotted green line) beyond D > 2, 
thereby confirming the regularizability of the ancilla dimension for a state of Gisin-Massar form. The missing data points in 
(a) and (b) reflect calculations beyond the lower bound of the numerical precision. 



III. REGULARIZATION OF THE ANCILLA 
DIMENSION 

The analytical study of Delgado et al already clarified 
that a multiqubit state of the Gisin-Massar form exhibits 
a rank deficient MPS representation with linear scaling of 
the bond dimension M with the number of clones. Nev- 
ertheless, the resource costs of the result yet seems unaf- 
fordable for large number of qubits. One may then won- 
der if any optimization protocol can be devised to allow 
the accomplishment of the same task with lesser amount 
of resources. In other words, we pose the following ques- 
tion: how well a state of Gisin-Massar form \GMm{^)) 
with bond dimension D (and equivalently requiring an 
ancilla of dimension D for its sequential generation) can 
be represented if only an ancilla with a smaller number of 
levels, D < D, is available? More formally: given a state 
\GMm(iP)) , with bond dimension D, what is the opti- 
mal MPS \GM M (tp)) of lower bond dimension D < D 
that "regularizes" the original huge ancilla demands by 
minimizing the distance 

_min || |GM M (V)) - |GMm(V)> f ■ (?) 

dim{\GM M W))=D<D 

We exploit two techniques that were developed by 
Saberi et al (2tj [38| in the context of sequential gen- 
eration of entangled states to perform the MPS approx- 
imation above, both exploiting a suitably designed local 
optimization of the ^-matrices in Eq. ([T]). In the first 
approach, use is made of a corollary of the singular value 
decomposition theorem from linear algebra to perform a 
local optimization procedure which might be referred to 
as "MPS compression", in analogy to the image compres- 



sion technique already used in computer science and engi- 
neering [39] . Let the SVD of matrix A with rank( A) = r 
be given by 

r 
i=l 

Then, the best possible lower-rank approximation to A 
that minimizes the Frobenius-norm distance (40| 

min \\A-A\\ F , (9) 

ranh(A)—r<r 

is given by [H, l35j 

r 

A = J2 a ^ v \ ■ ( 10 ) 

i=l 

The latter indicates a truncation scheme in which one 
retains only the f largest singular values of A to form the 
desired optimal lower-rank matrix A and discards the 
rest. Applying the outlined truncation scheme to each 
matrix VS, k = l,...,n, in Eq. <JTJ) , yields an MPS of 

lower bond dimension D = D — (r — f ) which is locally 
closest to the original MPS of bond dimension D in the 
Frobenius norm sense. It is noteworthy to mention that 
the precision of the method hinges to a great extent upon 
how well-decaying the singular values spectrum of the 
underlying matrices are. 

In the second approach (26l. [30L l38j. a DMRG-inspired 
variational optimization of ^-matrices is performed (2(1 
l3ll HH, l4H - i43l | to obtain the best possible approximation 
to \GMm ("0)) m the space of all MPS of lower bond di- 
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Table I. Results for regularization of the ancilla-qubit interactions 



number of qubits (n) 


3 


5 


7 


1 — T (without auxiliary unitaries) 


0.5000 


0.4730 


0.5684 


1 — T (with auxiliary unitaries) 


8.6938xl0~ 9 


3.7714xl0~ 4 


1.6600xl0 -2 



mension D < D of the form 

l 



\i n , ...,h) , (11) 



.ii=0 



by solving the minimization problem of Eq. ([7]) under 

the constant-norm constraint (GM m {4>)\GM m{iI))) = lj 
which is taken care of by introducing a Lagrange mul- 
tiplier A. Varying Eq. (|7|) with respect to the matrices 
defining \GM m(V0) leads to a set of equations, one for 
each ifc, of the form 



d 



(1 + A) {GM M W\GM M (ip)) 



-2Re{(GM M (^)|GM M (V'))} 







(12) 



which determines the optimal ^-matrices of the desired 
state \GM m (V 1 ))- These equations can be solved very 
efficiently using a "sweeping procedure" (30j in which one 
fixes all but the fc'th ^-matrix and solves the correspond- 
ing Eq. (fT2")l for the matrix VS. Then one moves on to 
the neighboring site and, in this fashion, sweeps back and 
forth through the chain until the convergence is reached. 
The cost function contains the overlap terms between 
the target state \GMm(iP)) and the simulational one 
GMm('0)} which can be straightforwardly calculated in 
MPS representation as illustrated in Fig. EJ^b). 

We have implemented numerically the two regulariza- 
tion methods outlined above and the results are illus- 
trated in Fig. [3J Focusing on the cloning of an input state 
of the form \ip) = -^(|0) + |1)), we define the fidelity of 
the procedure to be the overlap between the target Gisin- 
Massar state and the simulational one of the form 



F=(GM M {ip)\GM M (i>)) . 



(13) 



Strikingly, Fig.[3Ja) suggests that, for instance, a Gisin- 
Massar state of n = 13 qubits (associated with M = 7 
clones) can be simulated with almost perfect fidelity 
(1 — T « 10~ 16 ) with a regularized dimension of D = 3, 
and to a lesser degree (1 — T ks 10 -2 ) with D = 2. Note 
that, the analytical result of Delgado et al requires an 
ancilla of dimension D = 2M = 14 for realization of 
such a task, whereas our numerics allows the experimen- 
tal realization of SQCM with an ancilla of dimension only 



D = 3 which is considered to be amenable to the exist- 
ing physical setups. The performed procedure can be re- 
garded as a kind of "tomography" of Gisin-Massar state 
within the subspace of the constituent dimensions. As 
is evident from Fig. [21a), the dimension components D 
beyond D — - 3 make negligible contributions to the rep- 
resentation of the state of Gisin-Massar with n = 13. 

Figure [3{b) shows the scaling of the fidelity with the 
number of qubits n using an ancilla of dimension D = 2 
and D = 3. Since variational optimization approach 
allows for the feedback of information through several 
sweeps, it generally performs better than SVD trunca- 
tion scheme. The SVD truncation scheme is only a local 
optimization of information in the sense that each matrix 
is being optimized independently from the rest of the ma- 
trices, whereas the variational optimization, though be- 
ing itself a local optimization protocol too, makes up for 
such a locality by performing several sweeps and allowing 
the feedback of entanglement through the chain. 

All in all, the simulation promises the possibility of 
realizing SQCM with a manageable amount of ancillary 
levels Z? = 3upton = 15 qubits. We believe the reg- 
ularizability of the Gisin-Massar state mainly originates 
from symmetrization requirements of the optimal proto- 
col already discussed in previous section. 



IV. REGULARIZATION OF THE 
ANCILLA-QUBIT INTERACTIONS 

A state of Gisin-Massar form can be generated se- 
quentially provided that the required ancilla dimension 
D and ancilla-qubit unitaries are available pH [23l . l2t| . 
However, in general, some of the required local ancilla- 
qubit unitaries may not be accessible to a given physical 
setup. Given such a relevant constrained optimization 
problem fjij . we wonder if an optimization protocol can 
be devised by which a Gisin-Massar "target" state can be 
approximately generated with maximal fidelity. 

In general, the unitary time evolution of the joint sys- 
tem ancilla-qubit at step k of the sequential generation 
is described by a general unitary : Ha <8> T~Lq — > 
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U A ®U Q 



(14a) 
(14b) 



where Htjfi is a general bipartite Hamiltonian that en- 



tangles the ancilla with the fc'th qubit, hjj^j are real- 
valued coupling constants, and cti, <72, 03 are the usual 
Pauli sigma matrices, with <7q = I the identity matrix. 
For simplicity, we have considered the case D = 2, but 
similar generators of SU(iV) may readily be employed for 
D > 2. 

As an illustrative case, now suppose that only a re- 
stricted set of unitaries produced by the experimentally 
realizable ent ang ling Hamiltonian of the form of the 
XXZ-mode\ 45j are available 



#[fc] C = h[ i ] fa ® CT i + CT 2 ® era) + ftiffa ® ^3 , (15) 



containing two nonzero couplings h\ 



,[*] 



and 

Given a state of the Gisin-Massar form and 



'11 



the restricted entangling Hamiltonian of Eq. (|15[) . the 
aim is to tailor the restricted unitary operations fj£ft = 



'[fc] 



e " J i fc i in such a way that when applied sequentially 
to an arbitrary initial state of the joint system |$/) = 
\<Pi) ® yield a state of the form 



l*>=^i a -^ a ^] a l*/>. 



(16) 



that is "closest" to the target state of the form \(pp) <g> 
\GMm(iP)), where \<p F ) is arbitrary [4a |. In the ideal 
case, when the fidelity reaches unity, the ancilla can be 
set to decouple unitarily in the last step [23]. However, 
this ceases to be the case in general when the allowed 
ancilla-qubit unitaries are restricted. Thus, the optimiza- 
tion problem shall be carried out with respect to the cost 
function of the from 



/[!*)] =|| |<l) - \<p F ) ® \gMmW) \(' 



(17) 



involving a multivariable cost function in \<pf) and the in- 
teraction couplings {h^\ h 2 , ■ ■ ■ ; , h 2 ^}, as the vari- 
ational parameters, which can be solved in an iterative 
procedure developed by Saberi et al in the context of se- 
quential generation of entangled states [26J: One starts 
by picking a particular unitary, say U£P , and minimizes 

the cost function in Eq. (fTT)) . varying over {h^\ h^} and 
regarding couplings of all the other unitaries as fixed. 
Then one moves on to the neighboring unitary and op- 
timize its couplings. When all unitaries have been opti- 
mized locally, one sweeps back again and so forth until 
convergence. Varying over the vector \ip F ) an d using the 
resulting optimal one, the cost function boils down to 



/[|§)] = 2(1 -||(*|GM M (V))|| 



(18) 



(<Pf 



V 




V 

v [2] 












Figure 4. MPS graphical representation of the overlap pattern 
for optimization of ancilla-qubit interactions upon incorporat- 
ing the auxiliary single-qubit unitaries. 



suggesting the definition of the fidelity of the procedure 
as 



F=\\(V\GM M (i>)}\ 



(19) 



For the restricted entangling Hamiltonian of Eq. (fT5)) , 
however, the variational space is so small that the vari- 
ational optimization of the restricted unitaries per se 
does not result in much overlap with the target state 
\GMm(iP)), as illustrated in the first row of Table [I] 

However, we have found out that the fidelity can be 
improved upon systematically by allowing optimized ro- 
tations on the initial and final states of the qubits and 
ancilla, namely augmenting the restricted unitaries by 
unrestricted local "auxiliary" operations u^- and u A as 
depicted schematically in Fig. 01 In this way, our numer- 
ics imply that a faithful representation of Gisin-Massar is 
obtained with a restricted experimentally realizable en- 
tangling Hamiltonian of XX Z form up to n — 7 qubits, 
as the second row of Table U illustrates. 



V. CONCLUDING REMARKS 

In conclusion, we have found strong numerical evidence 
that a state of Gisin-Massar form can be faithfully rep- 
resented with an MPS of effective bond dimension D =3 
up to n = 15 qubits. Equivalently, it can be generated 
within sequential prescription with an ancillary system of 
experimentally manageable dimension D = 3. We have 
also analyzed the possibility of realizing SQCM with a 
restricted class of physical interactions and have demon- 
strated that such restrictions do not pose any obstacle to 
sequential realization of cloning machine. 

The current investigation and our results promises the 
possibility of elevating the SQCM ideas from an abstract 
theoretical level to practical recipes of experimental em- 
inence. They help to bridge the gap between theory 
and experiment for realization of a sequential "factory" of 
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quantum cloning under real-life restrictions and in a form 
amenable to the existing physical setups. The achieve- 
ments will be of importance, in particular, for realiza- 
tion of scalable quantum computing within all poten- 
tially scalable and hitherto realized QCMs, either in op- 
tical systems or NMR setups pWlEll. The method 
may accordingly prove useful for the recent experimental 
demonstration of probabilistic quantum cloning (47| in an 
NMR quantum computer. Generalization to general se- 
quential quantum cloning schemes for N — > M [48[ may 
be also the subject of further investigation. 
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